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Abstract 

Cochlea is an important auditory organ in the inner ear. In most mammals, it is coiled as a 
spiral. Whether this specific shape influences hearing is still an open problem. By employing 
a three dimensional fluid model of the cochlea with an idealized geometry, the influence of the 
spiral geometry of the cochlea is examined. We obtain solutions of the model through a conformal 
transformation in a long- wave approximation. Our results show that the net pressure acting on the 
basilar membrane is not uniform along its spanwise direction. Also, it is shown that the location of 
the maximum of the spanwise pressure difference in the axial direction has a mode dependence. In 
the simplest pattern, the present result is consistent with the previous theory based on the WKB- 
like approximation [D. Manoussaki, et al, Phys. Rev. Lett. 96, 088701 (2006)]. In this mode, the 
pressure difference in the spanwise direction is a monotonic function of the distance from the apex 
and the normal velocity across the channel width is zero. Thus in the lowest order approximation, 
we can neglect the existance of the Reissner's membrane in the upper channel. However, higher 
responsive modes show different behavior and, thus, the real maximum is expected to be located 
not exactly at the apex, but at a position determined by the spiral geometry of the cochlea and the 
width of the cochlear duct. In these modes, the spanwise normal velocities are not zero. Thus, it 
indicates that one should take into account of the detailed geometry of the cochlear duct for a more 
quantitative result. The present result clearly demonstrates that not only the spiral geometry, but 
also the geometry of the cochlear duct play decisive roles in distributing the wave energy. 
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I. INTRODUCTION 



The primary functions of the mammahan cochlea are to detect and analyze the acoustic 
signals in terms of their intensity and frequency contents. It is able to probe and analyze 
sound waves over wide ranges of frequency and intensity, e.g., for human being, over twelve 
orders in intensity and over three orders in frequency (20Hz - 20kHz) Some mammals 
are able to perceive frequencies as high as lOOkHz [2|. Remarkable capabilities appearing 
in hearing are essentially governed by both passive mechanical and active biophysical pro- 
cedures in the cochlea. The physics of hearing is pioneered by von Helmholtz in 1880s [s^. 
He proposed that different regions along the longitudinal length of the cochlea vibrate inde- 
pendently, working as a series of oscillators driven by the pressure of the fluid filled in the 
cavity of the cochlea. The coupling among the successive parts of the cochlea is thought to 
be mainly determined through the fluid, instead of the cochlea itself. Benefitted from the 
achievements in both the experimental techniques and the theoretical methods, auditory 
physics has been developing continuously from the time of von Helmholtz and expanding on 
ever- widening problems J Q, [s , Q] . 

Observations revealed that the cochleae are coiled in almost all mammals, monotremes 
being the exceptions. This auditory apparatus is a fully fluid-filled system of canals contained 
inside the petrous bone. It suggests that detailed calculation of the fluid motion and the 
pressure distribution within the cochlea should be crucial for understanding of hearing in 
any mechanical models of the cochlea. Fluid in the cochlea is a lymph. Physically, it can be 
treated as incompressible and inviscid as a reasonable approximation [?]. The cavity of the 
cochlea is partitioned by the basilar membrane (BM) into two channels which are connected 
through the helicotrema at the apex of the cochlea. At the other end, the base, the channels 
are sealed by two movable membranes working as pistons. These two movable membranes 
are labeled as the oval and the round windows, respectively. The channel abutting to the 
oval window is named as scala vestibuli (upper channel). Another channel abutting to the 
round window is known as scala tympani (lower channel). The upper channel is further 
separated by the Reissner's membrane (RM) into two subducts. However, the RM is usually 
ignored in the mechanical models for the cochlea because of its extreme flexibility, when 
the axial flow behavior is the main focus. As a consequence, the structure of the cochlea 
is typically modeled by two channels with the BM as their interface in most mechanical 
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studies. However, when the spanwise pressure distribution is included in the study, the RM 
affects the flow characteristics and, thus, cannot be neglected. We will address this problem 
below. 

The stapes is attached on the oval window and forms the boundary between the middle 
and the inner ear. Incoming sound waves are transmitted into the fluid in the channels of the 
cochlea through the vibration of the stapes. Nearly entire mechanical sound energy exerted 
by the stapes is used to produce the fluid movement. The sound pressure transmitted to 
the cochlear fluid interacts mechanically with the BM and vibrates it into oscillation in the 
form of a traveling wave that propagates along the BM from the base to the apex. The 
traveling wave on the BM grows as it travels, reaching a maximum at a position known 
as the characteristic place determined by the frequency of the stimulus, and then declines 
rapidly. The characteristic behavior along the longitudinal cochlea is described by a place- 
frequency map; high frequencies are resonant near the basal end, whereas low frequencies 
are resonant near the apical end jj]. 

The first direct observation of the responses of the BM to sounds through experiment 
was achieved by von Bekesey Q| . He measured the BM movement of the cadaver and found 
that the mechanical response is a traveling wave propagating from the base to the apex. 
Subsequent in-vivo measurements revealed that the responses of the BM are sharply tuned 
and very sensitive at low sound-pressure levels. Hence it is evident that hearing involves an 
active process. However, a passive mechanism also plays an essential role in hearing and 
also in understanding the auditory phenomena. 

Although the spiral coiling is one of the most distinct characters of the cochlea, it has 
been neglected in many mechanical models j^. Therefore, it is one of the prime interests 
to ask whether the shape of the cochlea influences hearing or just facilitates the packing of 
the long BM into a small space. If the answer is the latter, then the cochlea can be safely 
modeled as a straightly shaped duct to avoid the complexity of the spiral shape. Due to 
the mathematical difficulties in dealing with the shape of the cochlea, this problem has not 
been fully answered yet. Earlier studies mostly concluded that the shape of the cochlea has 



little effect on the response of the BM along the center line j^, lo, 10, 11|. However, recently, 
Manoussaki et. al. employed a WKB-like and linearly approximated solution to examine 
the effect of the curvature of the cochlea on the redistribution of the wave energy density. 
According to their linear expansion, it was concluded that the spiral geometry of the cochlea 
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does play an essential role in creating a tilt of the BM between the outer and the inner wall, 
which is believed to be important to produce the vibration of the hair cells. The calculated 
tilt is shown to be proportional to the curvature of the center line of the cochlea and, thus, 

n 

the maximum of the tilt is predicted to be always at the apex p^. In the present paper, we 
employ a three dimensional fluid model without the assumption of the WKB-like solution. 
By employing a conformal transformation and a long-wave approximation, we explicitly 
obtain the solutions of the fluid model which satisfies the boundary condition. 

Our analytic solution shows that the Manoussaki et al's solution 12j is naturally included 
as the lowest component of the general solutions, and has no flux across the channel width. 
The general solution shows that the resultant pressure acting on the BM is not uniform across 
the channel width. It is shown that different responsive modes show different patterns of 
the spanwise pressure tilt along the longitudinal length of the BM. Especially, for higher 
modes which allow the flux across the channel width, the magnitude of the pressure tilt is 
not always a monotonic function of the curvature as predicted by Manoussaki et al. jl2 |. 
Since the real fluid oscillation is a linear superposition of several modes, we expect that the 
the maximum of the resultant pressure tilt is not located exactly at the apex as predicted 



12j, but at a position determined by the input wave and the geometry 



by Manoussaki et al. 
of the cochlea. 

The arrangement of this paper is as follows. In Sec. II, we give the fluid model of the 
cochlea and calculate the velocity potential and the flow of the fluid in the duct of the 
cochlea. Sec. Ill discusses the pressure in the fluid using the Navier-Stokes equation. Sec. 
IV gives a discussion on the results. Sec. V concludes the paper. 



II. MECHANICAL MODEL 



Because most mammalian cochleae are coiled into a spiral, the cochlea is approximately 
modeled as a logarithmic spiral duct with a rectangular cross section and the RM is approx- 
imately considered as the center line in the upper channel as shown in Fig. 1, although the 
detailed geometry is much more complicated. We also assume that it stays in plane in the 
present model although the center line of the cochlea coils out-of-plane in reality. The BM 
partitions the duct into two channels, the upper and the lower channels. The height and 
the width of the cochlea are assumed uniform along its longitudinal length. In what follows. 
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we mainly discuss the response of the fluid in the outer part of the upper channel, between 
the RM and the outer wall using our idealized geometry. The result can be extended to 
the lower channel by conside ring the connectivity of the helicotrema and the symmetry as 



elaborated in the references 
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13|. Cochlea fluid is approximated as incompressible and 
irrotational, i.e., V ■ U = and V x U = 0, where U is the velocity of the fluid. Hence 
the velocity of the fluid in the cochlea can be described in terms of the velocity potential, 
U = V$. Thus, the velocity potential $ satisfies the Laplace equation 

V2<l> = 0. (1) 

In the cylindrical polar coordinate system, 

92 Id 1 92 92 

+ olJ7^ + 7^1^ + 1:^' (2) 



dR^ RdR R^dQ^ 

where R and G are the coordinates in the plane of the BM; C, is the coordinate perpendicular 
to the BM. At the apex, the polar angle is zero and at the spiral axis, i? = 0. Recently 
Manoussaki et al. studied the cochlea curvature effect on the radial shearing in the BM 
region 



12j. Instead of solving the Laplace equation, Eq. ([T]) directly, they assumed a WKB 
like propagating wave solution. 



$ = $o(i?,e,C)exp 



iujt-i I K{Q)Rmi<d)de 
Jo 



(3) 



where $0 is a slowly varying function of the angular coordinate B, Rm the distance to 
the spiral duct midline from the spiral axis, K the axial wave number, and u the input 
frequency. They applied this solution to the continuity and interface equations of the BM 
motion to obtain the BM displacement. In the first order approximation, they found that 
the displacement of the BM-fluid interface gets tilted across the width of the partition and 
the spanwise amplitude slope varies as l/Rm- Thus, they predicted that the maximum of 
the tilt to be always at the apex in the long-wavelength limit and the spiral shape of the 
cochlea plays an important role in distributing the wave energy density towards the cochlea's 
outer wall. 

However, it is believed that the pressure is transmitted through the fluid and not through 
the BM itself. Therefore, it is necessary to study the detailed fluid response to the applied 
pressure in order to obtain detailed information on the pressure redistribution and the 
sound propagation. In the present study, we solve the Laplace equation directly without the 



assumption of the WKB-like plane wave solution. Since the above equation does not contain 
any explicit time-dependence, solutions of this equation cannot describe propagation of the 
waves along the cochlea. However, it can give information on the spanwise distribution of 
the pressure on the BM of the steady response. It has been demonstrated that at least to 
the first order approximation, the in-phase assumption is a good approximation to discuss 
the spanwise profile of the BM Separating the variable, (, from the rest two variables, 
we have $(-R, 6, C,t) = (f){R, Q)H{()T(t). Eq. ([1]) is separated into two independent parts 

'^-^'Hia (4) 

dR^^ RdR^ R^dO'^ 

where /i is a separation constant that can be any constant. Here, we note that if we apply 
the WKB-like solution, Eq. ([3]), to Eq. ([5]), we obtain the same result as in Ref. Now, we 
impose a boundary condition that the wall of the cochlea is impermeable, i.e., the velocity 
of the fluid normal to the wall is zero on the walls of the cochlea. However, because the 
RM is extremely flexible, the velocity normal to the RM is not necessarily zero. Since the 
RM is fixed on the BM, the pressure distribution near the intersection line between the RM 
and the BM is expected to be much complicated. The solution for such a case can only be 
obtained through very complicated numerical calculations. Since, it is the aim of this paper 
to investigate the curvature effect of the cochlea and improve from the WKB-like solution, 
we neglect the edge effect and assume that the RM behaves as a free bondary. Therefore, 
the solution of the velocity potential in the vertical direction to the BM is 

H{0 = cosh [/i(C - h)] , (6) 



where h is the half height of the cochlea. For simplicity, we assume that /i is 



sufficiently 



small, /i <^ 1. Note that the small /i corresponds to a long- wave approximation [12, llSj- We 
will discuss later that the small /i is also related to a weak stimulation. Accordingly, Eq. ([5]) 
is approximately rewritten as 

dR^ RdR R^dQ^ ' ^ ^ 

The solution of Eq. ([7]) thus can be obtained as 

(/)(i?,e) = {Aicosie + BismlQ)R^, (8) 



where / is any integer and Ai{Bi) the constants to be determined by the boundary conditions. 

terms are neglected to avoid diverging contributions at small R. In order to employ 
techniques developed in the harmonic functions, the velocity potential, 0, is rewritten in a 
complex plane 

4>= [Ai{z' + z') - tBi{z' - z')] , (9) 

where z = Rexp{iQ) and z is the complex conjugate of z. The solution should satisfy the 
boundary condition, impermeability of the spiral walls. It thus requires to recombine these 
separated variables together in an appropriate way, so that the spiral boundary condition 
can be conveniently achieved. Conformal transformation is a natural way to achieve it. We 
introduce a complex potential, f = (p + iip. For / being analytic, and must satisfy the 
Riemann-Cauchy relation, 

^ dip 84^ ^ _dip_ . 
dx dy ' dy dx' 

where x and y are the real and the imaginary parts of the complex variable, z, i.e., z = x+iy. 
It is easy to verify that the imaginary part of the complex potential is the streamline of the 
fluid, ip = —[Ai{z'' + z}) + iBi{z^ — z^)]. The complex potential is thus 

f = z'[A-tB,]. (11) 

Now, we introduce a conformal transformation 

w = Poe-'\ (12) 

where po is chosen to be the distance from the origin of the spiral axis to the outer wall of 
the apex. After the conformal transformation, the complex potential has the form 

f=(^X^!^\ [A,-iBi). (13) 



Using that w = rexp[i6'] in the polar representation, 

I 



f 



-6 + iln — 

po 



(Ai-tBi). (14) 



The horizonal velocity potential 0(r, 6) is the real part of the complex potential Com- 
bining the horizonal and the vertical parts of the velocity potential, one has 

$(r, 9, C, t) = 0(r, 9)H{C)T{t) = Re{f)H{C)T{t), (15) 



where Re{f) is the real part of the complex potential /. Accordingly, the velocity potential is 
also a function of I. Different I corresponds to each different responsive mode. To distinguish 
different modes, we introduce a suffix I into the two-dimensional and the three dimensional 
velocity potential, velocity and pressure, i.e., (pi, Ui and Pi, respectively. In general. 



where [a;J represents the largest integer that is smaller than or equal to x. The above 
solution shows that the originally separated variables are recombined through the conformal 
transformation in the process of generating the spiral boundary. 

The position vector of the logarithmic spiral is r = pe^^ir in the polar coordinate system, 
where p is the distance from the spiral axis to the apical end of the spiral when the polar 
angle is zero; 6 is a parameter describing the spiral. The covering area of the spiral increases 
with b describing the compactness of the logarithmic spiral. On the outer wall, = poe^^Cr, 
the boundary condition of the impermeability of the cochlea wall requires that the normal 
component of the velocity should be zero; 



where e„ is the unit vector normal to the spiral (see Appendix A). This condition provides 

relations between Ai and Bi as follows for / = 1, 2, 3; 





(16) 



t/,,„(r„,^) = UKro,^)-e„ = 0, 



(17) 



Si = -hA 



(18) 



2b 



A 



(19) 



62 - 1 



■2, 



^3 = 



(3 - 62), 
362-1 



^3. 



(20) 



And the velocity potentials for the three lowest modes are 




(21) 




(22) 
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_ «J°(3,^_„,L)1. (23) 

1 - po V PoJ \ 

Since, we assumed that the RM is extremely flexible and, thus, behaves as a free boundary, 
the normal velocities of the inner side and the outer side are to be the same on the RM. We 
explicitly get the solution of the velocity potential of the outer part of the upper channel. 
According to Eq. fl2T|) . the normal velocity is always zero in the first mode. Consequently, 
this solution can be directly extended to the inner part of the cochlear duct. In this case, the 
fluid only has angular flux and the RM can be neglected in the mechanical model. However, 
for Eqs. (l22l) and ( |23ll . the fluid possesses normal flux. The normal velocity on the outer wall 
is zero but it is nonzero on the RM. Thus, it is expected that the RM critically influences 
fluid flow for the higher modes. In order to investigate this point more closely, we carried out 
a similar calculation on the inner part, from the inner wall to the RM. We find the normal 
velocities of the inner side and the outer side on the RM are virtually the same. Thus, we 
belive that when the stimulation is weak so that no excessive pressure can be build up at the 
RM and at the BM boundary, the assumption of the free boundary for the RM can provide 
at least a qualitative picture for the BM oscillation. 



III. PRESSURE IN THE FLUID ACTING ON THE BM 

The pressure of the fluid is governed by the Navier-Stokes (N-S) equation [17|. Viscosity 
of the lymph in the cochlear duct is usually negligible in the studies of hearing j^]. The 
inviscid N-S equation reads 

Pmass^ = -VP, (24) 

where pmass is the density of the fluid; D/Dt = d/dt + \J -V; P is the pressure of the fluid. In 
the study of hearing, the disturbance of the fluid in the cochlea is usually small. Therefore, 
the pressure in the fluid can be obtained through the linearized N-S equation [l7] 

The pressure in different modes can thus be obtained through the linearized N-S equation 
Pi = -ipmassUJ(pi{r, 9) cosh[/x(C - h)]T{t), (26) 
10 



where we assume an input frequency uj. The fluid pressure right above and below the BM 
are Pup and Pdown- When the upper channel is driven alone, the cochlear duct is expected to 
reduce the encountered large impedance through the lower channel. It has been shown that 
due to the symmetric structure of the upper and the lower channels and the connectivity of 
the helicotrema, we obtain an approximate relation P^p ^ — Pdown 0, 13] • Hence the net 
pressure in the mode /, which is the pressure difference between the upper and the lower 
channels acting on the BM, is 

Pf^{T,Q,t) = -2P,{T,Q,^,t), (27) 



where Pi is the pressure in the upper duct in the mode /. We deflne a dimensionless ti 
pressure on the BM, which can be directly compared to the BM displacement tilt in Ref. 

as 



t of 



12|, 



Pf^ {r,,e)-Pf ^{r,,e) 



^^l = oBM/^ o ^ ' K^^) 



where ro{6) is the distance to the outer wall following the spiral arc, 6'base = ^Ktt with K 
being the number of the spiral turns, and rc{6) the distance to the center of the BM. 

An arbitrary position vector on the BM can be expressed as r = + cTe„, where a is the 
normal distance from the outer wall to the position. Hence 

r=\r\= Tog i-^j , (30) 



\ ; o / \ 2 



where Ro{= VoVb"^ + 1) is the radius of the curvature of the outer wall. Then, the dimen- 
sionless tilt of the pressure can be rewritten as 

^pBM ^ Mro,0) - Mro9c,0) ^^2) 

where gc = 9 (^^) ^^"^ is the half width of the BM. For the flrst and the second modes, 
we have 



2 



(2K7r)2(62 + 1)2 



2 2 + 

In Qc ^^"^dc 



(34) 
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IV. RESULTS AND DISCUSSIONS 

Under the long- wave approximation, /i^ <^ 1 [12|, Il5|, we have obtained the solutions of 
the three dimensional mechanical model of the spiral shaped cochlea. Now, we show that 
fi can be sufficiently small for any frequency as long as the stimulation is weak. As an 
example, we study the case of mode = 2). Consider the velocity of the BM near C = 
at which the velocity of the BM reaches the maximum. The velocity of the BM is assumed 
to be the same as the velocity of the fluid on the BM. Therefore, the velocity of the BM can 
be represented by the velocity of the fluid as 

[/2,c(r, ^, 0, t) = -fj.^h^{r,e). (35) 



In a weak stimulation, the velocity of the BM is very small, ~ 10 ^m/s or less [18| in the 
passive system. Considering that the height of the cochlea is in the order of millimeter jl9|. 

fi (p[r,0) = ^ ~ 10 /s. (36) 

Hence, the small fi limit is consistent with not only the long-wave approximation but also 
the weak stimulation. 

In the case of the lowest mode (Z = 1) in which the flow has no normal flux, our results 
show that the dimensionless tilt of the pressure is a monotonic function of both the distance 
from the apex, 6, and the radius of the curvature, Ro- The maximum is at the apex, as 
shown in Fig [2l This is consistent with the previous theoretical result [12] . As shown in 
the Fig. [2)3, the dimensionless tilt of the pressure is a scaling function of the radius of the 
curvature, AP^^ ~ R~^'^^. This is in a close agreement to the result of the WKB-like 



approximation, AP^^ ~ [12|. In this mode, the solution can be directly extended from 
the outer part to the inner part. According to Eqs. (|2T|) and fl26|) . the pressure difference on 
the BM is also a monotonic function from the inner wall to the outer wall. In addition to this 
lowest mode, we find that there exist higher modes which show more complex behaviors. For 
example, in the second (/ = 2) and the third modes (Z = 3), the positions of the maximum of 
the tilt are determined by the compactness of the cochlea (described by b) and the reduced 
width of the cochlea {(Tc/ Po) as shown in Figs. [HHSl In these two cases, the solutions of the 
outer part cannot be directly extended to the inner part and the role of the RM may not 
be neglected. Our calculation reveals that except in the lowest mode, the distribution of 
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the resultant pressure acting on the BM is not only determined by the curvature of the 
cochlea but also influenced by the width of the cochlear duct (Fig. Hj), which again implies 
the importance of the RM geometry in a more quantitative calculation. The amplitude 
of the tilt of the first responsive mode is larger than the tilt of the second and the third 
responsive modes in the apical region (Fig. E]). 

The resultant pressure gradient in the radial direction induces the vibration of the hair 
cells on the BM. The pattern of the pressure tilt along the longitudinal length of the BM 
is therefore important for understanding the hearing mechanism. Previous work with the 
WKB-like approximation has shown that the curvature influences the tilt monotonically 12|. 
However, our results reveal that the influence of the shape can be more complex than that 
obtained by Manoussaki et al. [l^. 

The general solution of the pressure distribution in the fluid will be a linear superposition 
of various modes. To show the diversity of the pressure tilt along the BM, we now consider, 
for simplicity, only three lowest = (ai$i + a2^2 + «3'^'3)/(ai + 02 + cts), where ai, 02 
and 0:3 are any real constants. In order to determine the relation between different A;, we 
assume that the pressure at the basal end is same to the input pressure when each individual 
mode is allowed to exist alone. Since the coefficients appearing in the linear combination 
for a real oscillation cannot be determined a priori, thus, we assume two simple ratios of 
ai : 0:2 : 0:3 = 1 : 2 : 1 and ai : 0:2 : 0:3 = 2 : 2 : 1 and plot the tilt of the pressure of these two 
examples as superpositions of the modes in Fig. [7l The result shows that the maximum may 
not occur at the apex, but at a position not far from the apex. We believe that the exact 
maximum position is determined by the spiral and duct geometry of the cochlea and the 
detailed nature of the pressure wave. This result clearly indicates that the spiral curvature 
plays a decisive role in the wave energy redistribution in the cochlea as reported in the earlier 
calculation, although the detailed nature of the maximum position may be different from 
the WKB-like approximation Here, we note that we have assumed that the channels 
are fully filled with the fiuid. If the channels are partially filled, the formalae and te results 
will be totally different. 
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V. CONCLUSION 



In a long-wave approximation, the flow field of the three dimensional fluid is obtained 
analytically for the spiral shaped cochlear model. Pressure of the fluid is then calculated 
through the N-S equation. Result shows that the pressure acting on the BM is not uniform 
along the spanwise direction of the BM. Furthermore, the patterns of the pressure tilts 
along the longitudinal length of the BM are different for different modes. In general, the 
real pressure distribution is expected to be a linear superposition of various modes, although 
the simplest one may be dominant or at least substantial. Contributions from the higher 
modes are expected to push the pressure tilt maximum towards the basal end from the apex 
to a position which is determined by the curvature and the duct geometry. In conclusion, 
we have shown that the curvature of the cochlea plays a decisive role in hearing in addition 
to providing a compact storage. 
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APPENDIX A: UNIT VECTORS OF LOGARITHMIC SPIRAL 

Let the position vector of the curve is r. The unit tangent and normal vectors, as shown 
in Fig. [HI are respectively 



where is the unit vector in the r direction in the polar coordinate system. Consequently, 




(Al) 



In the polar coordinate system, the position vector of the logarithmic spiral is 



r = 




(A2) 



1 



{bir + eg) , 



(A3) 



e* = 



VWTT 



1 



{-ir + bee) ■ 



(A4) 



e. 



n 
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lower channel 



FIG. 1: Sketch of the shape of the cochlea in the present model. The top view of the cochlea is 
modeled as a logarithmic spiral; the cross section of the cochlea is approximated as a rectangle 
with constant width and height along the entire length of the cochlea. The duct of the cochlea is 
divided by the BM into two channels. 
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FIG. 2: (color online). The pressure tilt in the mode 1 {I = 1). The maximum always appears at 
the apical end. For ad Po^ we use 0.3. (a): The tilt is plotted as a function of the distance from the 
apex, (b): The dimensionless pressure tilt is a scaling function of the radius of the curvature, Rq- 
The slope of the straight lines is about 1.08. Therefore, it reveals that, in this mode, our result is 
consistent with the previous theoretical result [l^ . 
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FIG. 3: (color online). In the mode II (Z = 2), the tilt of the resultant pressure acting on the BM. 
It reaches a maximum and decreases along the longitudinal length of the cochlea. The position of 
the maximum depends on od Po and h. 
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FIG. 4: In the mode II, the maximum of the tilt appears at a position along the longitudinal length 
of the cochlea. It is determined by both the compactness and the width of the cochlea. 
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FIG. 5: (color online). The tilt of the pressure in the mode III (1 = 3). The curves are calculated 
at (Jc/ Po = 0.3. The maximum of the tilt appears at a position determined by the geometry of the 
cochlea. 
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FIG. 6: (color online). A comparison of the tilt of the resultant pressure for the first three modes. 
6 = 0.10, K = 3 and cJc/ Po = 0.3. In the apical region, the tilt of the first responsive mode is larger 
than the tilt of the second mode; The tilt of the second responsive mode is larger than the tilt of 
the third mode. 
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FIG. 7: (color online). Linear superpositions of the first three modes. The maximum of the tilt is 
located at a position determined by the coefficients of the superposed modes and the geometry of 
the cochlea, b = 0.10, K = 3 and dc/po = 0.30. Solid line corresponds to qi : 02 : as = 1 : 2 : 1. 
Dash line corresponds to qi : 02 : as = 2 : 2 : 1. 




FIG. 8: The unit tangent vector, it, and the unit normal vector, e„, of a logarithmic spiral. 



21 



